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Whereas single- and two-photon wave packets are usually treated as pure states, in practice they 
will be mixed. We study how entanglement created with mixed photon wave packets is degraded. We 
find in particular that the entanglement of a delocalized single-photon state of the electro-magnetic 
field is determined simply by its purity. We also discuss entanglement for two-photon mixed states, 
as well as the influence of a vacuum component. 



INTRODUCTION 



Consider an entangled state containing one or more 
photons. By how much is the entanglement degraded 
when the photon wave packets are described by mixed 
rather than pure states? For example, suppose Alice and 
Bob are given a two-photon polarization- entangled state 
[E 0, B A Hi — say the singlet state — but they are not 
told what the color of the photons is. All they know is 
the photons are either both blue or both green. Does 
their ignorance reduce the amount of entanglement they 
possess? The answer in this case is negative: the entan- 
glement is still one ebit, even though the overall state is 
mixed. After all, they could, in principle at least, apply a 
local measurement on each photon that measures its color 
but not its polarization. That way they are guaranteed 
to end up with a pure maximally entangled state. The 
fact that polarization and color are independent degrees 
of freedom is crucial here. 

Suppose now Alice and Bob are given a mode- entangled 
state jo, m M containing merely a single delocalized pho- 
ton falil M, G3, OS H 



beam splitter 



|0)a|1)s + |1)a|0)j 
V2 



(1) 



where A and B denote specific modes in Alice's and Bob's 
labs and |0) and |1) denote Fock states with zero and 
one photon, respectively. The notation used implies that 
modes A and B are well-defined. But suppose that Alice 
and Bob actually do not know what the color of the single 
delocalized photon is, green or blue (each with 50% prob- 
ability) . Or suppose they do not know the polarization of 
the photon, only that it is either left-hand or right-hand 
circularly polarized. Then they really should ascribe a 
mixed state to their field modes: an equal mixture of 
and the similar state 



|V/) := 



|0)A<|l)B' + |l)A'|0)i 

V2 



(2) 



where the primed modes refer to modes of different color 
or different polarization. In Section IIIII we will find the 
logarithmic negativity [HI, [l6[ of this mixed state to be 
Ejsf = log 2 (1 + \/l/2) < 1, so that in this case Alice's and 
Bob's state does lose some of its glamour [a pure state 
of the form |T]) would, obviously, contain one ebit of en- 
tanglement]. Note that, indeed, Alice and Bob cannot 




(|1>|0> + |0>|1» 



FIG. 1: A single photon impinges on a 50/50 beam splitter 
and one ebit of entanglement is created between the two out- 
put ports. How does the entanglement change when the input 
photon is in a mixed state? 

Answer: When we use the logarithmic negativity [lEl Il6| 
as our entanglement monotone of choice we find £V = 
log 2 (l + ^purity) , in terms of the purity Trp 2 of the input 
state (see Section iHl]) . 



use the same local filtering measurement of frequency to 
filter their state to a pure entangled state: as soon as a 
photon is detected on, say, Bob's side, the state collapses 
to either |0)a|1)b or |0)a'|1)-B'i an d not to the desired 
pure entangled state or \ip'}. Alternatively, a nonlocal 
filtering measurement of color could upgrade Alice's and 
Bob's state to a pure entangled state, but that nonlocal 
operation could and would increase the amount of en- 
tanglement. The distinguishing feature of this example 
compared to the previous one is that color or polariza- 
tion and mode are dependent degrees of freedom. More 
precisely, color and polarization are part of what defines 
a mode. (Note, by the way, that we can, of course, cal- 
culate the logarithmic negativity of the state in the first 
example as well. The result (see Appendix) is that it 
indeed equals unity.) 

The purpose of this paper is to continue investigat- 
ing questions of this sort: by how much is the entan- 
glement of single- or two-photon states degraded when 
the photon wave packets are not pure but mixed? See 
Figs. [T]JJ]for typical examples of questions considered in 
the present paper. The motivation for this research is, 
of course, the simple fact that typically any photon pro- 
duced in an experiment is represented by a mixed state 



2 



beam splitter 



|1) 



k 



%/2 



|2)|0)-|0)|2)) 



FIG. 2: Two identical photons impinge on a 50/50 beam split- 
ter: what is the entanglement at the output? 
Answer: 1 ebit, similarly to the answer in Fig[T| the output 
is a delocalized photon pair. 



beam splitter 




(|1>|0> + |0)|1)) 



([1'>|0)-|0)|1'» 



FIG. 3: Two distinguishable photons (with, say, orthogonal 
polarizations or different colors) impinge on a 50/50 beam 
splitter: what is the entanglement at the output? 
Answer: 2 ebits (2 equivalent versions of an entangled delo- 
calized single-photon state) . 



@, [13, 0, EH- For instance, even if one's source pro- 
duces a Fourier-limited wave packet, in practice one will 
not know exactly the timing of the wave packet, or the 
exact central frequency, or the exact width. For another 
example, consider a single photon heralded by detection 
of the other photon of a down-converted pair of photons. 
Whenever there is some entanglement between the two 
photons, tracing out one photon necessarily leaves the 
remaining photon in a mixed state. Our study comple- 
ments those in Refs. [2(| where the effects of spectral 
entanglement on polarization degrees of freedom are in- 
vestigated. 

In Section |TT] we start out by collecting some useful 
results about the description of single- and two-photon 
wave packets to be used in later sections. In Section UTTl 
we focus our attention on entangled states that can be 
generated by splitting a single photon on a 50/50 beam 
splitter. Ideally that leads to an output state with one 
ebit of entanglement, but, as we will show, the entan- 



FIG. 4: Two distinguishable photons with orthogonal po- 
larizations impinge on a 50/50 beam splitter. When there 
is already entanglement between the polarization degrees of 
freedom in the input state, what is the entanglement at the 
output? 

Answer: .Eout = 2 + Ei n /2 for the entropy of entanglement, 
generalizing the answer illustrated in Fig. [3] (see Section ITV)) . 
How does this answer change when the input state is mixed? 
Answer: it's complicated . . .see Section ITVl for details. 



glement resulting from splitting a mixed single-photon 
input state is less than one ebit but turns out to be a 
simple function of its purity. In the same Section we will 
also consider the more realistic case of a nonzero vacuum 
component of the state of the field and its effect on entan- 
glement. In Section [IV] we consider typical mixed states 
of two orthogonally polarized photons arising from type- 
II down conversion, and in that case too we consider the 
effects of the presence of a (large) vacuum component. 
We also compare the mode entanglement one obtains by 
splitting both photons on a 50/50 beam splitter with the 
entanglement that may already be present in the input 
state between two orthogonally polarized modes, and find 
one may increase the amount of entanglement that way. 
We conclude the paper with a summary that also dis- 
cusses some possible extensions of the present work. 



II. SINGLE- AND TWO-PHOTON WAVE 
PACKETS: PRELIMINARIES 

A. Single photons 

Consider a single photon of a definite polarization 
propagating in a well-defined direction. Then a pure 
state can be described in terms of continuous modes [22| 
as 



1, 



dtip(t)a\t)\v) 



(3) 



where |v) is the vacuum state, a*(t) an operator that cre- 
ates a photon at time t and ip(t) the temporal mode func- 
tion of the wave packet. Often it is more useful to Fourier 
transform this representation into frequency space, and 



3 



describe the same state by 



B. Two photons 



1^,) = / dujtjj(uj)a'(u})\v). 



(4) 



Typically, the function ^(w) will be appreciable only in 
a small bandwidth a around a central frequency u>q, with 
it « wo, and the wave packet ([4]) describes a quasi- 
monochromatic photon. Since the creation operators 
bear the relation 



t(t) 



1 



dusa/ (. 



to e 



(5) 



the spectral shape ip(u)) is thus completely determined 
by the temporal mode ip(t) in the sense that 



1 



(6) 



The single-photon states displayed so far are pure. Mixed 
states of single photons arise, e.g., when they are part of 
a multipartite system and we trace over the other par- 
ties; or if one of the transverse degrees of freedom of the 
photon is traced out; or if we are simply ignorant about 
one or more of the properties of the photon. In any case 
the mixed state of a single photon is described in terms 
of the density matrix 



Pi 



dXP(X)\l^)(l^ 



(7) 



Here A stands for any parameter or combination of pa- 
rameters that can possibly be involved in the mode func- 
tion, although no concrete form is given yet. Some good 
examples of what A could stand for are arrival time 
("time jitter"), central frequency ("frequency jitter"), the 
width of the (Fourier-limited) wave packet, etcetera p^j . 
When writing the mixed state in terms of parameters 
with such a clear physical meaning, the states |l^, A )'s 
will in general not be orthogonal. However, we can al- 
ways diagonalize the density matrix, and rewrite it as a 
discrete sum involving orthogonal states, 



Pi = X^ fe l lfe )( lfe l 



Here p k and |lfc) are eigensolutions to 



with 



pi | lfc) = Pk\lk) 



Ojk, 



(8) 



(9) 



(10) 



and ^2 k Pk — 1- For our purpose of calculating entangle- 
ment of single-photon states, the latter representation is 
often more useful. 



Now consider states of exactly two photons. Since we 
have in mind the two-photon component of the state pro- 
duced by type-II down conversion [lj], [H, HiJ , we assume 
the photons have orthogonal polarizations. We will sim- 
ply indicate the ordinary and extraordinary polarizations 
by "H" and "V" . Pure and mixed states for such photon 
pairs can then be expressed as 

|2 V; ) = I duj [ d^w,w>HM°v( w/ )|v>. (11) 



and 



P2 = / d\P(A)|2^><2^|, 



(12) 



respectively. The commutators of the mode operators are 

[a k {u),a\,{<J)]=6 k , v 8(u>-wr), (13) 

where {k, k'} = {H, V}. 

It is well-known that for a pure state of a bipartite 
system there is always a discrete Schmidt decomposition. 
In our case this allows us to rewrite 



|2v>) = £Va^UI 



(14) 



For a state like Eq. ([lip we can explicitly define the new 
creation operators as 



,t 



wit = / dw</>fe(w)a| A (w), 



(15) 



where A&, ip k an d <fe can be obtained by solving the eigen- 
value problems [23, 



duj' pa{uj, cj')cpk(u>') = Afc(^fc(w), 
dijj'pB{u,u')4>k{u') = \k<l>k{u), (16) 



with the "reduced density matrices" given by 

pa(u J ,u i ) = J i/^(m/)f( u V'), 

Pb{u,u') = duj"ifj(ui" ,lu)i(;*(uj" ,u>'). (17) 



The mode operators h\.'s and v^s satisfy the standard 
commutation relations 



(18) 
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C. Mode entanglement 

In the rest of the paper we will calculate mode entan- 
glement between field modes 0, H| that are spatially 
separated [27|. For that purpose we expand the den- 
sity matrix in an appropriate orthonormal basis of Fock 
states, including states with no photons, single-photon 
states |lfe) as defined above, etc. 

We use two standard measures of entanglement in 
this paper: one useful entanglement measure, but only 
valid for pure bipartite states, is the entropy of entangle- 
ment [28[. For example the mode entanglement between 
modes of orthogonal polarization of the pure state \2^) 
of Eq. (jlip is defined in terms of the Schmidt coefficients 
appearing in Eq. (|14[) as 



(19) 



E{\2 4 ,)) = -^A fc log 2 A fc . 



The other entanglement monotone we will make use of 
is the logarithmic negativity (l5l [l6j , defined in terms of 
the partial transpose (PT) of a matrix. For example, for 
a two-photon system PT is easily defined if we expand its 
density matrix in the basis spanned by |fc)i|Z)2i('^|2(^|, 
with \k)i, \l}2, \tti)i and \n)2 referring to orthogonalizcd 
single-photon states respectively. Then we have 



pt r 
P -> P , 



r 

Pklmn 



Pknml i 



(20) 
(21) 



Pklmn being the matrix element. We then solve for the 
eigenvalues of p r , ending up with a series of real e^'s, 
since both p and p T are Hcrmitian. The logarithmic neg- 
ativity [H, HU is defined as 



Ejs{p) = log 2 ||p r ||i = log 2 



ek\ 



(22) 



It is worthwhile to notice that the absolute sum of all 
eigenvalues J2k \ e k \ equals 1 plus twice the absolute value 
of the sum of all negative eigenvalues of p v . 



III. MODE ENTANGLEMENT OF 
SINGLE-PHOTON STATES 

A. No vacuum component 

1. General results 

When an optical field containing exactly one photon in 
a pure state is split on a 50/50 beam splitter, the output 
state looks like 



1 

72 



-~(|0>|1> + |1>|0» 



(23) 



and possesses exactly one ebit of entanglement [12, Il3l - 
[ljj]. If, on the other hand, the photon entering one in- 
put port of the beam splitter is mixed, the calculation of 



entanglement in the output state is a little more compli- 
cated, but the problem can still be solved analytically. 
We consider single polarization only and start by an in- 
put state that is already expanded in its diagonal form, 
Eq. © 



5>*IW*I- 



(24) 



fc=i 



We put it on a 50/50 beam splitter, which in the Heisen- 
berg picture transforms the mode operators for the input 
ports a* , w into those for the output c 1 , $ in the follow- 
ing way 



1 

71 



i i 

i -i 



(25) 



In this picture the state stays unchanged, but when ex- 
pressed in terms of creation operators at the output ports 
it looks different: 



Pout 



Ph 



^Pfc|lfc)a ®a (lfc| 
fc=l 

1 " 

-^p fc (|l fe ) c |0) d + |0) c |l 



k/d) 



k=l 

®(c<lfc|ii< 
n 

k=l 



(26) 
(27) 



where each submatrix pk can be written out as 



Pk 



/0 

Pk /2 Pk /2 
pk/2 p k /2 

\o o oo 



(28) 



Here rows and columns correspond to states |0) c |0)d, 
\0) c \Wd, |lfe) c |0)d, |lfe)c|lfe)d and their conjugates re- 
spectively. Naively, taking the PT simply gives 



Pi 



p k /2 

pk/2 

pk/2 

Pk/2 



(29) 



In the total density matrix context, the two pk/2 on 
the diagonal remain independent of other /5fc's. These 
give rise to nonnegative eigenvalues of p T . Care must 
be taken, however, with the two off-diagonal elements. 
Those matrix elements share the vacuum state |0) c |0)d 
with all other submatrices, and therefore beg an eigen- 
value solution to the matrix 



pi/2 
Pi/2 

. Pn/2 



Pn/2 






(30) 
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It can be shown, with a little effort, that the only two 
nonzero eigenvalues of this matrix are 



A1.2 - ±- 



Notice from Eq. © that 



£p8 

k=l 



1 



Erf- 



\ fe=i 



Trpl 



in' 



(31) 



(32) 



and that Trp 2 n is an invariant quantity under basis trans- 
formations, or equivalently, unitary operations. A^a can 
thus be evaluated as 



\ 1 , 2 = ± 1 -(Trp? n )^. 



(33) 



Since we learn from the definition that ||/9 r ||i equals half 
of the absolute value of the sum of all negative eigenval- 
ues, we immediately get 



EAf( Po nt) = log 2 [l + (Trp: 



2 sl/21 



(34) 



as announced in Fig.[TJ (Of course, the purity of the out- 
put state equals that of the input state, as we assumed 
the beam splitter to be describable by a unitary opera- 
tion.) 



2. Example 

With this conclusion, we attempt to calculate the log- 
arithmic negativity for single-photon mixed states for 
which all wave packet functions ip(u>) are Gaussian- 
shaped 



, (w-w ) 2 
ip{u!) oc exp [ 



(35) 



and which is mixed with respect to the arrival time r 
of wave-packet peaks. Assuming a Gaussian distribution 
for arrival times as well, we write 



P(t) cx exp ( - 



So, explicitly the density matrix is 

2 



(36) 
(37) 



Pin = A dr exp 



<7- 



du> / du>' 



exp 



(W - Ll>o) 2 (w' - UJof 



xa'(u)|v) (v|a(a/), 



+ i(u' — ujq)t 

(38) 



with a normalization coefficient A. Since the state is 
profiled in spectral space, we have 



Trftn = / duj(v\a(uj)p in a\u))\v} 



(39) 



By using the commutation relation 

[aH.aV)] =5(w-a/) 

and basic algebra we can show that 

Trpl = {l + Aa^l)-^. 



(40) 



(41) 



The logarithmic negativity after the beam splitter is then 
according to ([34]) 



EjsfiPoMt) = log 2 [l + (1 + 4<tV; 



(42) 



which in cither limit a r — or a = reduces to unity. 
The first limit is just a pure state, which is easily con- 
ceived, while the latter fact is a bit harder to reveal, but 
think of ip{oj) as S(ui — u> ) when a — > 0. Then, according 
to Eq. (|38p but with a different normalization constant 
A', we have 



= A' dr exp [ - 



du>8(u) — u) )e 



dw'S(uj' — w )e luJ T a) (lo)\v) (v\a(iv') 



pm 



= A' J drexp f-_j oT(w )|v)(v|o(w ), (43) 

which is equivalently a pure state. Physically, S(u> — uj ) 
corresponds to monochromatic light whose wave function 
extends homogeneously along the time axis to both in- 
finities, so that the concept of wave-packet arrival time 
no longer applies. 

In conclusion, the entanglement of the output state 
depends only on the ratio of the "incoherent" width of 
the mixture in time, <j t to the "coherent" width in time 
of each wave packet \/a. For a large incoherent width 
the entanglement reduces to zero, as expected. 



B. Adding the vacuum 

Real experiments involving single photons typically are 
described by a state involving a vacuum component in 
addition to the single-photon component. The phase be- 
tween the vacuum and a particular single-photon Fock 
state may or may not be known or controlled. We may 
write a state containing the two Fock states as 



sCVT^vl + vfcj-^il), (44) 

with p the fixed a priori probability for the state to con- 
tain exactly one photon. We consider two extreme cases 
here: 

• f(f) ~ when the state is pure and we simply 
represent it, instead of p, as a pure state 



|vaci in ) = V 1 - pW) + VpI 1 ); 



(45) 
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• f(<p) is a flat distribution so that the cross terms 
vanish after the integration. The state is thus re- 
duced to 

Pvacl ln = (l-J>)|v)(v|+p|l)(l|. (46) 

We treat these two cases one by one. The output state 
of Eq. (f4"5")) after a 50/50 beam splitter is 



|vaclout) - VT^\00) + J|(|10) + |01)). (47) 

The entropy of entanglement and the logarithmic nega- 
tivity are straightforwardly calculated and the result is 

£(|vacl out )) 
= i-I[(i + v /i- p 2 )log2 ( 1 + A / 1 _ p 2) 

+ (1 - V^P 1 ) log 2 (l - y/T=?)\, (48) 

^(|vacl out )) =log 2 (l+p). (49) 

The latter expression is particularly simple. Of course, 
both measures of entanglement vary between and 1 for 
p varying between and 1. It can be shown that Efj is 
always larger that E when p e [0, 1]. 

For the mixture, we are at liberty to assume the diag- 
onal expansion Eq. ((8]) along with Eq. |2| and Eq. (|10|). 
Hence we have 

Pvacl out = (1- P )|00)(00| 

+f E A *(l 01 *> + M»« 01 fcl + (lfc0|)(50) 
fc 

The quantities we are interested in are now 

JMpvacw) = (l-p) 2 +p 2 J2^l (51) 

k 

Eu(p™x J = \og 2 { P + Pur 1 ' 2 ), (52) 

where Pur denotes the purity of the input state. We see 
that Eq. (|52p generalizes the expression Eq. (|49)l which 
is only a special case when Pur = 1. Of course, it also 
generalizes Eq. (|34[) . In conclusion, E^j- depends only on 
p and the purity, which in turn is also affected by the 
value of p. 

IV. ENTANGLEMENT OF TWO-PHOTON 
STATES 

A. No vacuum component 

A broadband-pumped down-conversion process pro- 
duces a state whose two-photon component can be writ- 
ten in the form pip , with the mode function in the ideal 
(pure-state) case described as [17( 

ip(uj ,uj e ) — a(uo + w e )$(w„,w e ). (53) 



Here w Gje are the frequencies of the two photons with 
ordinary and extraordinary polarization, respectively. 
a(u) +u} e ) is the pump spectrum envelope, and the phase- 
matching function $(w , to e ), after a great deal of simpli- 
fication, is 

3>(w -I- V ,U) e + v e ) 
= svac{[v {k' -k' p ) + v e {k' e -k' p )]L}, (54) 

where v a%e are deviations from the perfect match fre- 
quencies dj ,e, and where L is the length of the nonlinear 
medium. Moreover, k' p e are the first derivatives of wave 
vectors with respect to frequency for pump photon and 
outcoming o- and e- photons respectively at the perfect 
phase-matching condition 

uj p — 0j o + Q e . (55) 

We still restrict ourselves to Gaussian wave packets only, 
and so we assume 

a[y Q + v e ) ~ cxp I — I , (56) 

so that ip(cj ,u) e ) is real everywhere. In general, due to 
our ignorance about the precise timing of the pump pulse 
or about its precise central frequency, the state generated 
will actually be a mixed state. 

We are interested in the entanglement of such a (pure 
or mixed) state between the two orthogonally polarized 
modes. Moreover, just like in the preceding section, 
we wish to calculate the entanglement that results from 
splitting such a two-photon state on a 50/50 beam split- 
ter. The resulting entanglement after the beam splitter 
is of a different sort, it's entanglement between the two 
output modes of the beam splitter, not between orthog- 
onal polarizations. An interesting question is whether 
that entanglement is larger or smaller than the initial 
polarization entanglement. 

Although we will have to resort to numerical meth- 
ods to calculate both types of entanglement, we can an- 
alytically determine the relation between pure-state en- 
tanglement before the beam splitter and that after the 
beam splitter: Assume that the Schmidt decomposition 
of a state [those coefficients can be obtained, in some ap- 
proximation, analytically [29| . but we won't need them 
explicitly] described by Eq. ([53]) is 

12m) =53^4^1^, (57) 
k 

(Here we have associated o-photons with horizontal po- 
larization and e-photons with vertical polarization.) The 
entropy of entanglement and logarithmic negativity are 
thus expressed in terms of Schmidt coefficients as[15j 

E(|2 in )) = -^A fe log 2 A fe , (58) 

k 

S Ar (|2 in )) = 21og 2 (E\/A^)- (59) 

k 
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It can be shown that the state after the beam splitter 
can be Schmidt-decomposed similarly as 

|2 out ) = i|2) c |0) rf -i|0) c |2) d 



E 



fx 



*i.t „,t 



L ck u dk\ 



E 



■«^S*|v>, (60) 



- after BS 
before BS 



where h' ck , v' ck , hl k and trL are associated with c^ H {ui) 1 
Cy(uj), d\j(uj) and d v (uj) according to Eqs[T5j \2)c,d 
stands for a specific two photon state, take c e.g., 



|2)c = E^L«Liv) 

which in frequency space is actually 



(61) 



|2) c = / / dv e a(y + i/ e )$(i/ , i/ e ) 

4(^)4(^e)|v). (62) 

The notation for the state is just a shorthand notation 
emphasizing its orthogonality with respect to any single- 
photon state or vacuum state. Hence we may conclude 

1 



£(|2 out )) = 2 + -£(|2 in )), (63) 
£AT(|2 out )) = 21og 2 (l + 2 E ^l 2 -))/ 2 ). (64) 



One obvious question is now whether these same relations 
will still hold for mixed input and output states. 

Since analytical solutions to the two-photon problem 
seem impossible without further approximations, even for 
the pure-state case (but for an exception see Ref. j2?|), 
we therefore turn to numerical methods where we intro- 
duce some standard approximations that were used be- 
fore in [25l [26| . Integrals over continuous frequency are 
converted into sums over discrete frequencies, and infin- 
ity as the integral limit is replaced by an artificial cutoff 
according to 



/OO />00 n n 
duj I duj' -> AwAw' ^2 E 
-OO J —OO • 1 1 1 



(65) 



3=1 k=i 



For convenience we hence choose Ato — Alu' and n = n' . 
The choice of Aw is determined by requiring the integrals 
to converge numerically. Moreover, we have to choose a 
scale for the many frequencies that occur in this problem. 
Quite arbitrarily, we have chosen to rescale all quantities 
with a dimension of frequency to Q, defined through the 
relation 



We also use 



(k' - k' )m = 2.25. 



(k' e - k' )m = 0.63. 



FIG. 5: Entropy of entanglement of pure two-photon states 
described by Eqs. (|53p - (|56jl before and after a 50/50 beam 
splitter as a function of the dimensionless pump width a/Q. 
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These two relations are similar to those used in Ref. [2f 



FIG. 6: Logarithmic negativity of pure two-photon states de- 
scribed by Eqs. (|53|l - (|56[) before and after a 50/50 beam split- 
ter as a function of the dimensionless pump width a/Q. 



In Figs. [5] and [5] we plot the numerically evaluated en- 
tanglement for pure states, as a function of the width a 
of the Gaussian pump pulse in units of f2. We verified the 
validity of Eq. (|6"3")l and Eq. (j6"4"| . Both v Q and v e (devia- 
tion from the perfect-match frequencies) are cut off from 
— 20 to 2D, . This is not driven by questions of numeri- 
cal convergence, but by the freedom one has to consider 
only those photons in a certain frequency interval. For 
instance, if one uses narrow-band detectors then only the 
entanglement between photons of frequency within that 
bandwidth will be relevant. We simply used cut-off val- 
ues such that the central peak and two side peaks of the 
sine function (|54[) are taken into account. 

Next we consider mixed two-photon states in a way 
that is similar to what we did for single photon states. 
We assume Gaussian wave packets with an uncertainty 
in arrival time. That is, we introduce a time-displaced 
mode function for two frequencies by 

i/)(w, J) -> ^ r (w, J) = Lu')e- tiuJ+u ' )T , (66) 

and we choose P(X) in Eq. (jX3J) to take exactly the same 
form as Eq. (|3"T|) only that A is replaced by r. Since the 
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FIG. 7: Logarithmic negativity of mixed two-photon states 
(|67p - (|68p after a 50/50 beam splitter vs. the dimensionless 
mixed-state width a T Q. The dotted curves are calculated by 
plugging the input state's Ej\f into Eq. (|64p and the solid 
curves are direct numerical results. Each set of dotted and 
solid curves has different wave packet width a. From top to 
bottom: a/Q = 0.5, 1,2. 



1.5 



0.5 











2p after BS 






2p before BS 






1 p after BS 















0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

Purity 



FIG. 8: Logarithmic negativities of various mixed output 
states of a 50/50 beam splitter are plotted against purity. 
The single-photon result follows directly from Eq. (f34)) , while 
the two-photon case is calculated from state Eq. (|67p . We 
have chosen a = fl for both two-photon states here and in all 
remaining figures. 



3.4 r 



Gaussian shape for r falls off fairly quickly, we decided 
to extend the integral in our numerical calculations over 
the interval (— 2cr T , 2a T ). To be more explicit, after the 
50/50 beam splitter the mixed state looks like 



P2out 



drP(r)|2 out (r))(2 out I 'J | •) 



(67) 



where P(r) ~ exp(— t 2 /<t 2 ) and 



1 



|2 ou t(V)) = - / dv a I dv e a T {v + z/ e )$(i/ , v e ) 



,-i(v +v<.)r / t 



( c LK) c v(^) + 4,(i/ e )4rOo) 
-c] I (u )d\ r (i' e ) - d]j(i/o)d\r(v e ))\v). 



(68) 



By realizing that the two-photon states can be expanded 
in their own space \2{t)) c ^ instead of single-photon com- 
bination space, the size of the density matrix is greatly 
reduced from (N + l) 4 x (N + l) 4 to (2N + T + l) 2 x 
(2N+T+1) 2 , N and T being the number of discretization 
intervals along u and r axes respectively. 

Fig. [7] shows our numerical results regarding mixed 
two-photon states. We display two kinds of curves for 
three sets of parameters. The solid curves give the nu- 
merical results for E^/(p out ), whereas the dotted curves 
plot 21og 2 (l + 2 B - V te»>)/ 2 ). In the case of a pure state 
these two quantities would be the same, according to 
but for mixed states there is a small difference, indi- 
cating the relation (|64p is a fair approximation for mixed 
states. 

From the preceding section we know the precise rela- 
tion between entanglement generated by a beam split- 
ter and purity for single-photon states. For comparison 
we plot in Fig. [5] the logarithmic negativity as a func- 
tion of the purity of the mixed state for both single- and 
two- photon states. The almost linear behavior indicates 



3.2 - 
3 - 

2.8 

2.6 

2.4 

2.2 
2 - 



from EQ 
- right after BS 
■ filtered after BS 



FIG. 9: Logarithmic negativities for the output of mixed two- 
photon states (|67p - (|68p . The dotted line is calculated from 
the relation Eq. (|64p . which is true for a pure state only. The 
dash-dotted line is the entanglement that would result after 
local filtering of the output state (see Section [IV CP . 



the close relation between the characteristic purity of the 
system and the amount of entanglement that can be ex- 
tracted under ideal conditions. Obviously, entanglement 
increases with purity. 



B. Adding the vacuum 

Type-II down conversion does not produce a two- 
photon state. Instead it produces a superposition of a 
two-photon state and the vacuum. Thus, let us first con- 
sider a pure state of the form 



|vac2 in ) = y/l-p\y) + VpJ2 V^*4 w fcl v )' ( 69 ) 

k 

where p is the a priori probability to detect two photons, 
which typically will be small. In terms of the Schmidt co- 
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efficients we can easily calculate the entanglement present 
in the state (jB"9"j) . 



-(1-p) log 2 (l -p) -p\og 2 p 
~P^A fc log 2 A fc , (70) 



£(|vac2 in )) 



3v(|vac2in» = 21og 2 ( v /T _ ^+ V / A^)-(71) 



Just as before we can express the entanglement of a state 
resulting from splitting this state on a 50/50 beam split- 
ter in terms of the entanglement of the input state. The 
50/50 beam splitter accordingly converts the state into 

|vac2 out ) = v/l - ^|0)| c 0) d + ^(|2) c |0) d + |0) c |2) d ) 



^EVm4/4+>4A)|v>, (72) 



where |2) is given by Eq. (fBTjl . The entanglement mea- 
sures are calculated to be 

J B(|vac2 out )) 



l-p/2 + VT^ l-p/2 + VT^ 
n lo S2 o 



2 l0g2 2 

+P - f lo g 2 P ^ f X] Afe l0g 2 Afe ' 



(73) 



^(|vac2 out )) = 21og 2 (l + VpEv / ^)- (74) 

fe 

We observe that 



E ut — 



l-p/2 + VT^P, l-p/2 + VT^ 

o lQ g2 X 



o 10 S2 o 



and 



1-p fl-p 

lo S2 



(75) 



(76) 



These relations are the generalizations of Eq. (|63|) and 
Eq. (|64|) respectively when vacuum is involved. 

For mixed states involving the vacuum no analytical 
results seem to be possible, so we reverted to numeri- 
cal calculations. Fig. [TO] plots some results from those 
calculations as a function of the probability p. 




FIG. 10: Logarithmic negativities for the output of a mixed 
two-photon state accompanied by vacuum at a fixed pump 
width a = and a mixture width a T Q, — 1. p is the propor- 
tion of the two-photon state, as indicated in Eq. (|69p . The 
dotted line is calculated from the relation Eq. (|74p . which 
is true for a pure state only. The dash-dotted line is the 
entanglement after the beam splitter after the local filtering 
operation discussed in the text (Section IIVC[) . 



C. Local filtering 

The mixed two-photon state (|67p is written as a mix- 
ture of pure two-photon states (|68p . The latter states 
are superpositions of two types of states: those with an 
odd number of photons in each output port of the beam 
splitter (namely, one), and those with an even number 
(namely, zero or two). By imagining performing a quan- 
tum nondemolition measurement of the parity of the pho- 
ton number on (one of the) output ports, we are applying 
a local filter. This filtering cannot increase the entan- 
glement and thus the entanglement after filtering gives 
a lower bound to the total entanglement present in the 
mixed state (IrJT)) . In the pure-state case we can calculate 
the amount of entanglement resulting from filtering by 
collapsing the output Eq. (f?2"]) into either even-number- 
photon state 

, . . 1 



(VT=p|0)| c 0)«j 
■^|2) c |0) d + ^|0) c |2) d ) 



(77) 



with probability 1 — p/2 or the odd-numbcr-photon state 

l^dd) = ^ £ v%(vl k hlk + hl k v f dk )W) (78) 

with probability p/2. Averaging the entanglement over 
the two possible measurement outcomes yields 



E 



(even/odd) 



= (l-p/2)E M {\cj) evei 
= P + plog 2 E 



))+p/2E N {\^ oAd )) 

iM'-i: 



1 



.(79) 
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Figs. [9] and [10] illustrate this filtering effect for mixed 
states as a function of er T f2 for fixed a — Q, and as a 
function of p for fixed a and cr r , respectively. One sees 
that about a third of the entanglement in the output 
state arises from coherence between states like |1, 1) and 
1 , 2) + 1 2, 0), in symbolic notation. 



V. SUMMARY 

We have quantified the entanglement of various mixed 
states containing exactly one or exactly two photons, 
as well as states with a nonzero vacuum component, 
by using the logarithmic negativity. For pure states 
we found simple relations between the entanglement of 
single-photon and two-photon states before and after a 
50/50 beam splitter. For mixed states such relations are 
still found to be approximately true. 

The simplest result arises for a mixed delocalized sin- 
gle photon. Its entanglement depends only on its purity. 
The result illustrates that even a perfect deterministic 
single-photon source (never producing more than a single 
photon) still may not be sufficient for certain quantum- 
information processing purposes (quantum computing 
based on dual-rail encoding, for example) if a large degree 
of entanglement is needed. 

We considered fairly realistic cases by explicitly includ- 
ing a vacuum component of photon states, as well as 
including the spectral and/or temporal shapes of pho- 
ton wave packets. But an obvious generalization of the 
present work would be to include full three-dimensional 
mode structures. Moreover, in the case of two pho- 
tons impinging on a beam splitter we only treated the 
case of photons with orthogonal polarizations (having in 
mind type-II down conversion), but the similar case of 
identically-polarized photons is interesting as well, and, 
perhaps surprisingly, more complicated. 



— horizontal(H) or vertical(V) — , but equally and clas- 
sically mixed in color — green(G) or blue(B) — , can be 
expressed in modes as 



1 



1 



where 



P = M(<f>i\ + ^\h)(M, 



i) = —(\GH) a \GV)b - \GV)a\GH) b ), 



and 



\4>2 



1 

V5' 



\BH) a \BV)b-\BV)a\BH)b). 



Expanding p in the basis of (BH, BV) A ® (BH, BV) B 
(GH, GV) A ® [GH, GV) B gives 
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The partial transpose of p in the same basis is then 
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Appendix 

Here we calculate explicitly the logarithmic negativ- 
ity for the first example to show the entanglement of a 
color-mixed polarization-entangled state is still one ebit, 
in spite of the mixed nature of the state. The two- 
photon singlet state, maximally entangled in polarization 



The eight eigenvalues are easily found to be this: six 
times the eigenvalue 4, and twice —\. This yields the 
logarithmic negativity of the state: 



£AT = log 2 (6x J+2x I- ~|) 



as we announced. 
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